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1  General  Introduction 

It  is  well  known  that  3D  woven  composites  have 
excellent  damage  resistance  and  damage  tolerance 
due  to  the  presence  of  z-yarns  [1].  According  to 
previous  experimental  researchers  3D  woven 
composites  shows  better  damage  tolerance  than 
laminated  textile  composites  without  z-yarns  such  as 
plain  woven  composites  even  though  there  is 
stiffness  and  strength  loss  due  to  z-yarns  [2]. 

Analytical  methods  and  numerical  methods  have 
been  widely  used  to  evaluate  the  effect  of  z-yarns. 

Analytical  methods  based  on  beam/plate  models 
predicted  the  apparent  or  effective  fracture 
toughness  in  the  presence  of  translaminar 
reinforcements  [3].  Finite  element  method  using 
discrete  spring  elements  and  cohesive  elements 
successfully  simulated  the  damage  behavior  of 
transversely  reinforced  composites. 

However,  recent  experimental  studies  have 
indicated  that  transverse  cracks  developed  inside  fill 
tows  (90  degree  tows)  induce  interlaminar 
delamination.  In  fact  this  phenomenon  is  very 
common  in  crossply  laminated  composites.  This 
characteristic  was  found  in  both  plain  woven  and  3D 
woven  laminates  at  early  stage  of  loading.  Hence,  it 
is  necessary  to  understand  and  characterize  the 
interaction  between  transverse  cracks  and 
interlaminar  delamination  before  investigating  the 
effect  of  z-yarn. 

In  this  paper  we  present  FE  modeling  of  quasi¬ 
static  short  beam  shear  test  of  plain  woven  laminated 
composites  (see  Fig.  1).  Cohesive  elements  were 
used  in  regions  where  transverse  cracks  and  multiple 
delaminations  were  expected  to  occur  based  on 
experimental  observations.  Parametric  studies  with 


various  properties  of  the  cohesive  elements  were 
conducted.  Predictions  of  peak  load  and  load- 
deflection  curves  including  post  damage  regime 
obtained  through  the  FE  simulation  were  compared 
with  representative  experimental  results. 

2  Modeling  multiple  delaminations 

Short  beam  shear  testing  set  up  and  FE  model  are 
shown  in  Fig.  1(a)  and  (b),  respectively.  Note  that 
only  half-length  was  modeled  because  of  symmetry 
with  respect  to  yz  plane.  Each  of  the  layers  was 
homogenized  as  an  orthotropic  material.  For  the 
purpose  of  reducing  computational  time  plane  strain 
elements  were  used  for  short  beam  shear  test 


(a)Setup  for  short  beam  shear  test 


y  ® - >X 

(b)FE  model  for  short  beam  shear  test 
Figure  1.  Experiment  setup  and  FE  model 


The  specimen  used  for  this  study  was  a  plain  woven 
laminated  composite  stacked  in 


(0^745^)2/0^ 


sequence.  Superscript  of  pw 


denotes  plain  woven.  Each  of  layers  consists  of  S-2 
glass  fibers  and  SC- 15  epoxy  matrix.  Material 
properties  for  0  degree  plain  woven  were  obtained 
from  the  previous  research  of  Xiao  et  al  [4].  These 
properties  are  shown  in  table  1. 


Table  1.  Elastic  properties  for  plain  woven 
composites 


Ei=E2 

E3 

G13=G23 

G12  V  13=V23 

Vl2 

(Mpa) 

(Mpa) 

(Mpa) 

(Mpa) 

27.5 

11.8 

2.14 

2.9  0.4 

0.11 

The  properties  for  45 -degree  plain  woven  layer 
could  be  obtained  using  appropriate  coordinate 
transformation.  In  order  to  model  intralaminar 
failure  (transverse  cracks)  and  interlaminar  failure 
(delamination)  cohesive  elements  associated  with 
damage  functions  were  employed.  The  cohesive 
elements  were  placed  between  layers  and  the  center 
of  fill  tows  for  interlaminar  delamination  and 
transverce  cracks,  respectively.  These  locations  were 
based  on  experimental  observations  where  most 
transverse  cracks  appeared  within  fill  tows  (90 
degree  fiber  direction).  Since  warp  tows  and  fill 
tows  comprising  woven  composites  are 
unidirectional  fiber  composites,  this  observation 
does  not  violate  a  general  rule  that  transverse  tensile 
strength  of  unidirectional  composites  is  less  than  that 
of  matrix.  And  it  is  reasonable  to  assume  that 
transverse  cracks  form  at  the  center  of  fill  tows 
since  specific  location  inside  the  fill  tows  have  no 
significant  difference  in  global  deformation.  It  was 
assumed  that  direction  of  pre-defined  crack  paths  at 
top  and  bottom  layers  were  in  the  transverse 
direction  (perpendicular  to  the  longitudinal  direction 
of  the  beam)  while  crack  paths  at  other  layers  were 
at  45 -degree  to  the  longitudinal  direction.  This  was 
based  on  the  principal  stress  direction  as  depicted  in 
Fig.  1(b). 

The  use  of  cohesive  element  associated  of  damage 
function  makes  it  possible  to  demonstrate 
progressive  delamination.  One  of  the  key  issues  is 


define  the  parameter  for  damage  function  of 
cohesive  elements.  Traction  separation  law  has  been 
used  for  damage  function  for  cohesive  element.  Bi¬ 
linear  cohesive  model  was  adopted  for  current  study. 
This  function  enables  stiffness  loss  to  be  calculated 
as  delamination  propagates. 


cj  =  {\-D)Kd 


(1) 


0  ,d<d0 


D  =  \ 


Q 

J?  ^  Q 

1 


,d0  <d  <df 
,df  <d 


where  a  is  traction,  K  is  stiffness,  D  is  damage 
variable,  d  is  displacement,  d0  is  displacement  at 
damage  initiation  and  df  is  final  displacement . 
Quadratic  stress  based  failure  criterion  was  used  for 
damage  initiation  for  both  transverse  crack  and 
interlaminar  delamination. 


where  <j°n  and  cj°s  are  interfacial  normal  and  shear 
strength. 

a 


Figure  2.  Bi-linear  damage  curve  for  cohesive  element 


Table  2.  Properties  for  cohesive  elements 


K 

< 

< 

GIC 

G„c 

(N/mm3) 

(Mpa) 

(Mpa) 

(N/mm) 

(N/mm) 

Transverse 

crack 

106 

50 

70 

1 

3.8 

Interlaminar 

delamination 

106 

30 

40 

1 

3.8 

In  addition,  quadratic  relationship  was  used  for 
damage  propagation  under  mixed  mode  of  normal 
and  shear  directions. 
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The  results  from  FE  model  are  shown  in  Figs.  3  and 
4.  From  Fig.  3  one  can  note  that  the  presence  of 
transverse  cracks  ensures  delamination  initiation 
consistent  with  the  experimental  observations.  Load- 
deflection  curve  shown  in  Fig.  4  is  helpful  in 
understanding  how  local  damage  patterns  affect  the 
global  deformation  behavior  of  SBS  specimens. 


Figure  3.  Comparison  of  damage  patterns 


Figure  4.  Comparison  of  load-deflection  curve 


delamination  occurs,  load  begins  to  decrease.  With 
more  deflection,  transverse  crack  and  interlaminar 
delamination  continue  to  develop. 

3  Parametric  studies 

3.1  The  effect  of  Mode  II  fracture  toughness 

Parametric  studies  were  performed  for  the  sake  of 
investigating  the  effect  of  inplane  strength  and 
interlaminar  strength  and  fracture  toughness  on 
global  behavior.  By  varying  the  cohesive  parameters 
for  transverse  crack  and  interlaminar  delamination,  it 
is  possible  to  demonstrate  the  effect  of  variation  of 
strength  and  fracture  toughness.  To  begin  with,  the 
variation  of  Mode  II  interlaminar  fracture  toughness 
values  (3.4  N/mm,  3.8  N/mm  and  4.2  N/mm)  was 
considered  since  Mode  II  fracture  is  dominant  under 
short  beam  shear  loading.  Material  properties  other 
than  Mode  II  interlaminar  fracture  toughness  were 
kept  constant.  Figure  5  shows  the  results  for  various 
Mode  II  fracture  toughness  values.  It  is  seen  that 
Mode  II  fracture  toughness  plays  an  important  role 
in  the  post  damage  regime. 


180 


— Experiment 
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Figure  5.  Comparison  of  load-deflection  curve  for 
various  interlaminar  fracture  toughnesses 

Lines  from  the  origin  to  the  point  corresponding  to  a 
deflection  of  2  mm  were  plotted  to  estimate  the 
residual  stiffness  of  specimens.  Higher  values  of 
Mode  II  fracture  toughness  lead  to  higher  residual 
stiffness  and  hence  higher  damage  tolerance. 


Transverse  crack  due  to  tensile  stress  at  bottom  layer 
at  early  stage  of  loading  cause  stiffness  loss.  Other 
transverse  cracks  developed  in  middle  layers  result 
in  interlaminar  delamination.  Once  interlaminar 


3.2  The  effect  of  interlaminar  shear  strength 

Interlaminar  shear  strength  is  another  important 
parameter.  Three  values  of  interlaminar  shear 
strength, 38,  40  and  42  MPa,  were  considered  in  this 


3 


study.  The  results  for  various  strengths  are  shown  in 
Fig.  6.  As  expected,  higher  value  of  interlaminar 
shear  strength  resulted  in  higher  peak  load  and  lower 
residual  stiffness. 


Figure  6.  Comparison  of  load-deflection  curve  for 
various  interlaminar  shear  strengths 

3.3  The  effect  of  transverse  crack  strength 

The  global  behavior  of  plain  woven  laminated 
composite  is  affected  by  interlaminar  delaminations 
which  are  induced  by  transverse  cracks.  So  it  is 
necessary  to  investigate  the  variation  of  transverse 
crack  strength.  In  this  study,  normal  and  shear 
strengths  were  changed  in  the  same  ratio  with 
respect  to  the  original  values.  Referring  to  Fig.  7, 
stiffness  drops  are  delayed  and  residual  stiffness 
tends  to  increase  as  transverse  crack  strengths 
increases.  However,  peak  loads  are  still  controlled 
by  interlaminar  shear  strength.  Strong  tows  prevent 
intralaminar  cracks  and  delay  delamination  initiation 
and  thus  can  increase  damage  tolerance. 


4  Future  work 

Three-dimensional  woven  composites  can  be  used 
for  enhanced  damage  tolerance.  In  order  to 
understand  the  effect  of  z-yarn  on  delamination, 
simulation  of  short  beam  shear  tests  of  3D  woven 
composites  with  10%  z-yarn  will  be  performed.  In 
the  simulation,  the  properties  for  cohesive  elements 
estimated  from  baseline  specimen  and  discrete  beam 
or  spring  element  to  model  z-yarns  will  be  used.  The 
results  from  simulation  will  be  analyzed  in 
conjunction  with  experimental  results. 


. '  Fill  tow 
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i 
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Fig. 2.  10  %  z-yarn  woven  composite  [2] 

5.  Summary 

This  paper  focused  on  the  simulation  of  short  beam 
shear  tests  on  plain  woven  laminates.  FE  model 
which  includes  transverse  cracks  and  interlaminar 
delamination  shows  good  agreement  with 
experimental  observations.  Parametric  studies 
showed  high  values  of  interlaminar  shear  strength 
caused  high  peak  loads,  and  increase  of  Mode  II 
interlaminar  fracture  toughness  resulted  in  increased 
damage  tolerance. 


Figure  7.  Comparison  of  load-deflection  curve  for 
various  transverse  crack  strengths 
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